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We propose a single crossover phase diagram applicable to 2D coUinear Heisenberg antiferromagnets 
(AFMs) and ferromagnets (FMs), and show that the scaling regimes of AFMs and FMs are in one- 
to-one correspondence. The phase diagram is split into classical and quantum regions. Our two key 
results are: (i) in the classical region, the AFMs and FMs exhibit nearly identical behavior near 
their respective ordering wavevectors, which we observe for S=l and higher using series expansions; 
and (ii) in the quantum region, quantum critical (QC) regime is present not only for the AFMs, but 
for FMs as well. 



I. INTRODUCTION 

We report a study of two-dimensional, coUinear, spin-S 
antiferromagnets (AFMs) and ferromagnets (FMs) aimed 
at better understanding the magnetic behavior of both, 
and in particular the role of quantum versus classical fluc- 
tuations in these systems. A number of quasi-2I? exper- 
imental systems, including spin-i AFMs La2Cu04 and 
SraCuOaCla 0gl , spin-1 AFMs La2Ni04 and K2NiF4 
II, spin- 1 AFM Rb2MnF4 |§ as well as spin-i FMs 
such as K2CUF4 §1 are well described over a range of 
temperatures by the 2D Heisenberg model 



(1) 



on a square lattice, which we study by high temperature 
series expansion methods. We expect much of our results 
to apply to other lattices having coUinear long-range or- 
der at T = 0. 

Let us begin by considering the relevant energy scales 
in the problem. One important energy scale for both 
AFMs (J > 0) and FMs (J < 0) is the T = spin 
stiffness ps, defined as rigidity with respect to a twist in 
the magnetic structure. Quantum 1/S corrections make 
this quantity different for AFMs and FMs with the same 
value of spin; for model (^, the order of magnitude is 
however the same, p^^^^ - JS'^, p™ = JS'^ (note that 
the FM value is exact). 

As shown by Chakravarty, Halperin, and Nelson ||] 
for AFMs, and by Kopietz and Chakravarty for FMs, 
the asymptotic T — s- magnetic behavior for these mod- 
els obeys scaling and depends on only two dimensionful 
parameters: the energy scale ps, and a quantity which 
sets the overall length scale and can be obtained from 
the g — > limit of the spin wave dispersion £{q). For 
AFMs, e{q) is linear and one defines the T = spin wave 



velocity as c = limq^o e{q) / q ^ JSa. For FMs, is 
quadratic and one defines the T = spin wave stiffness 
= limq^o e{q) / q'^ ^ JSa^. Note that the spin wave 
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stiffness p has the dimension of energy x length , while 
the spin stiffness ps has the dimension of energy. 

For the Heisenberg model, the spin wave spectrum is 
well-defined for all q and its upper bound, reached at the 
boundary of the Brillouin zone where q ~ 1/a, can be 
estimated as A ~ c/a ~ JS for the AFMs, or A ~ p/a^ ^ 
JS for the FMs, i.e. not only ps, but also A ~ JS has 
the same order of magnitude in AFMs and FMs. 

Whereas the T ^ behavior depends on ps and c for 
AFMs, or Ps and p for FMs, the T — 00 behavior (the 
Curie- Weiss law) depends on pcL = JS{S + 1) and the 
lattice constant a. In what follows, we show that this sim- 
ple observation is part of a general picture where the low 
temperature behavior, which may contain several scaling 
regimes, at T '--^ A crosses over to the high-temperature 
behavior, which may also contain several scaling regimes. 
A similar crossover in 3D occurs below the long range or- 
dering temperature (which is zero in 2D), and was studied 
by Vaks, Larkin, and Pikin 1^. 

In model (|l|), the ratio A/ps ^ 1/S can be made ar- 
bitrarily small by increasing S, but it cannot be made 
arbitrarily large as > ^- Nevertheless, the case 
A/ps 3> 1 is of great interest. First, there exist mod- 
els where ps goes to zero while A does not, in which 
case A/ps 3> 1 applies rigorously; an example of such 
a system is the two-layer Heisenberg model near 
the critical point where the Neel AFM long-range or- 
der vanishes. Second, for the S = 1/2 AFM model (Q), 
the ratio A/ps ~ 10 ^ 1 is quite large; furthermore, 
there exists evidence JTT| , p^p [ p^JTo| ] , first pointed out by 
Chubukov and Sachdev, that at intermediate tempera- 
tures this model is indeed in the quantum critical limit 
A/Ps > 1- 
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FIG. 1. A phase diagram of spin-S quantum antifer- 
romagnets (AFMs) and ferromagnets (FMs); see Table I 
for the corresponding scahng and crossover expressions for 
the correlation length. All regime boundaries are grad- 
ual crossovers rather than phase transitions, and their po- 
sitions can only be defined within numerical factors of or- 
der unity. The behavior in the T/A ^ 1 region is in the 
universality class of 2D classical magnets, where magnetic 
properties of AFMs and FMs for the same value of S are 
the same near their respective ordering wavevectors, and 
depend on T/pcL, where pcL = JS{S + 1). The classi- 
cal behavior includes Curie- Weiss regime (T S> max(ps,A)) 
and classical scaling regime (A <g T ^ Pcl), which are 
separated by a fairly wide classical crossover regime for 
A ^ T ~ pCL, where scaling does not hold, but never- 
theless ^AFM ~ Cfm ~ aipcL{T/JS{S + 1)). In the region 
T/A <C 1 quantum effects are important and therefore AFMs 
and FMs behave difi'erently. For the AFMs, this region is in 
the universality class of the QNLcr model. Here these scal- 
ing regimes are found for both AFMs and FMs: renormalized 
classical (RC) regime for T <^ min(ps,A); quantum critical 
(QC) regime for ps ^ T <C A; and a wide quantum crossover 
regime for T ~ pa <§; A, where the correlation length re- 
mains a universal function (different for AFMs and FMs) of 
the thermal de Broglie wavelength for spin waves (At ~ c/T 
or At = {p/Tf/'^) and the ratio T/ps. 



Regime 


Correlation Length 


Curie- Weiss 

T > max(pcL, A) 


Cafm — Cfm ^ a 
Properties depend on T/pcL 


Classical Crossover 
T - PCL > A 


^AFM — CfM 

= a?/'CL(T//3CL) 


Classical 
PCL > T > A 


Cafm = ^FM a{T/pcL) 
X exp(27rpcL/T) 


Quantum Critical 
Ps < T < A 


^AFM ^ c/T, 

^FM-ip/Tl0g{T/ps)f^ 


Quantum Crossover 
T ~ p, « A 


S.AFM = (c/T) (j)AFM{T/ps), 

Cfm = {p/ry/^MT/ps) 


Renormalized Classical 
T < min(ps, A) 


Cafm ^ {c/ Ps) exp(27rps/T), 
^FM - (pr/p,)i/2exp(2^p3/T) 



Table I. Correlation length in the scaling and crossover 
regimes shown in the phase diagram of Fig.|^, after Refs. 
[p^p|j7|jll|] , and this work. The sign ~ indicates the presence 
of a numerical prefactor. The functions (?!>afm and i^fm are 
universal; the function i/^cl is not, for instance, it explicitly 
depends on the ratio of the next-nearest and nearest-neighbor 
exchange couplings. 
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II. CURIE- WEISS REGIME 

At high temperatures T 3> max(ps , A) , the spins are 
weakly correlated and one can keep only a few lead- 
ing terms in the high temperature series expansion. 
For model (0), this temperature range corresponds to 
T 3> max( JS", JS*^). The corresponding mean-field the- 
ory yields the familiar Curie- Weiss law. 



III. RENORMALIZED CLASSICAL REGIME 

In the opposite hmit of T <C min(ps, A), the behavior is 
of the renormalized classical (RC) type. This regime was 
studied in detail in Refs. P,[7|,pT[ , where magnetic prop- 
erties were calculated by mapping the system onto the 
classical nonlinear cr-model and setting the momentum- 
integration cutoff to be proportional to the thermal de 
Broglie wave vector qx, defined such that e{qT) ^ T. The 
predicted ^(T) has the same form for AFMs (J > 0) and 
FMs (J < 0), when expressed in terms of their respective 
thermal de Broglie wavelengths At = ^/qt- 



e T ( 2'Kps 



(2) 



(here the exact value of the prefactor is obtained from the 
correlation length of the classical 0(3) nonlinear-cr model 
in the minimal subtraction scheme [TslJT^ ). However, not 
only the values of ps differ for FMs and AFMs, but also 
At has qualitatively different temperature dependences: 



At — 



c/T for AFMs 
y/]^ for FMs 



(3) 



An important question is upto what temperatures 
should this RC expression hold? We recall that when 
T is larger than the upper bound of the T = spin 
wave spectrum, A, there can be no spin waves with wave- 
lengths that are comparable to the thermal de Broglie 
wavelength. Since such spin waves are important for 
RC theory, one would expect it to be applicable only 
for T <C min(A,ps), which for model (|l]) translates into 
T ~ min( JS", JS^) (here, we ignore numerical factors of 
order unity). At higher temperatures new physics must 
arise. 

The Curie- Weiss regime occurs for T ^ max(ps , A) 
and the RC regime for T ^ min(ps,A). These two 
regimes exist for all S. We now turn to the discussion of 
the regimes where the temperature is larger than one of 
the scales ps or A, but smaller than the other. 
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FIG. 2. A semi-log plot of C/a vs. r/( J5'(5'+ 1)) for AFMs 
(solid symbols) and FMs (open symbols) for 5 = 1/2 (circles) 
and S — 5/2 (diamonds). We also plot the data from Ref. 

for the classical 5 — » oo model (solid line). For S = 1/2, 
AFM and FM correlation lengths are different and neither 
agrees with the S ^ oo limit. For S = 5/2, AFM and FM 
correlation lengths agree both with each other and with the 
S oo model. The data for S = 5/2 is plotted as represen- 
tative for large spin models; for all studied spins 5 > 1, the 
behavior is similar although the agreement between ^afm and 
^FM for S = 3/2 and S = 2 is not as striking as for S — 5/2. 
This data collapse provides evidence that any quantum ef- 
fects, which are expected to be different for AFMs and FMs, 
become unimportant already for fairly small spins 5 > 1, in 
agreement with the conjecture due to Ref. |17|. 



IV. CLASSICAL AND CLASSICAL CROSSOVER 
REGIMES 

In our earlier work in collaboration with Grevcn and 
Birgeneau (l^, we proposed a scaling crossover sce- 
nario to describe substantial deviations from RC behav- 
ior observed in the neutron scattering experiments for 
AFMs with spin-one and larger. This scenario calls for 
a crossover to the behavior of the 5*^00 classical sys- 
tem when spin wave energies for all wavevectors in the 
Brillouin zone become smaller than the temperature, i.e. 
for T > A. The correlation length at this RC to classical 
boundary, obtained from Eqs.(^), is exponentially large 
for large spin. 



exp(5') 
S 



> 1 for S* > 1. 



(4) 



In ||T^ , the arguments in favor of RC to classical crossover 
were based on data collapse for large spin when plotted 
versus T/{JS{S + 1)). The studied values of spin were 
small enough {S = 1/2 to S' = 5/2) to make the results 
sensitive to the choice of the renormalized temperature 
as T/{JS{S+1)) or T/JS^. 

Here, we resolve this arbitrariness by studying 
antiferro- and ferromagnets simultaneously. AFMs and 
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FMs have the same classical limit, but the quantum ef- 
fects differ and therefore the difference between AFMs 
and FMs can be used as a probe for the strength of 
quantum effects. It turns out that for any S" > 1, in 
the temperature range of interest ~ 10^), the correla- 
tion length calculated by high temperature expansions is 
nearly the same for AFMs and FMs with the same value 
of S (Fig.|^). Therefore, it must be determined by clas- 
sical rather than quantum physics. We emphasize that 
such an argument does not rely on any particular way of 
collapsing the data, because it is drawn from comparisons 
of AFM and FM models for the same spin. On the other 
hand, for S = 1/2 the difference between ^afm and ^fm 
is always large and temperature-dependent, which indi- 
cates the importance of quantum physics for this value 
of S. Spin-one appears to be a borderline case. 

Having established that in the temperature range T 3> 
A finite-spin models behave similarly to the S ^ oo clas- 
sical magnet, we now comment on the behavior of the 
latter. The asymptotic low-temperature behavior of the 
classical 2D Heisenberg antiferromagnet is given by : 

(5) 



a 



exp , 

32 27r/9cL V 



T 



where pcL = JS{S + 1). It turns out that because of 
the very small numerical value of the prefactor (~ 0.01), 
this formula becomes accurate only for very low tem- 
peratures, T/pcL ^ 0.6 > 100). We label the 
intermediate temperature range where neither Curie- 
Weiss nor classical scaling behavior apply the classical 
crossover regime. In this regime, neither large nor small- 
temperature asymptotic expressions describe the correla- 
tions accurately, nevertheless, all properties of the model 
are nearly the same for AFMs and FMs, and depend on 
T through T/pch only. 



V. QUANTUM CRITICAL AND QUANTUM 
CROSSOVER REGIMES 

We call quantum critical that behavior which can for- 
mally be obtained for ps <C T ^ A. For the antiferro- 
magnets, this regime is the asymptotic quantum critical 
(QC) regime The correlation length ^ is propor- 

tional to At = c/T with corrections that depend uni- 
versally on T/ps, and to leading order the dominant fre- 
quency scale for spin fluctuations is (D ~ T. 

For the FMs the spin wave spectrum is quadratic, 
which causes infrared {q — > 0) log divergences for ps <^ 
T ^ A. Such divergences (not present for AFMs) are 
cutoff by Ps , and lead to the following multiplicative cor- 
rection to the correlation length: 



FM 



{p/Tf'\og-'r^{T/p,), 



1/2 



(6) 



in either the sigma model or the Schwinger boson formal- 
ism. 



Note that such a behavior cannot be the true quantum 
critical behavior because the limit — > leads to singu- 
lar ^ at finite T. Therefore, Eq.(H) must fail as ps 0, 
and the T = quantum ferromagnet-paramagnet phase 
transition is likely to belong to a different universality 
class [|l^. Nevertheless, lacking a better name and for 
compatibility with the AFM nomenclature, we call this 
behavior a ferromagnetic quantum critical regime. 

The applicability of the QC description to model (|l|) 
for S* = 1/2 AFM is widely discussed in the literature 
1^,^ 12 ^,|l^. A similar analysis for the ferromagnets 
requires detailed calculations of their properties in the 
QC regime, and is beyond the framework of this paper. 



VI. CONCLUSION 

In this paper, we present a complete phase diagram 
of spin-S antifcrromagnets (AFMs) and ferromagnets 
(FMs). We show that comparing AFMs and FMs al- 
lows one to elucidate crossover effects from the quan- 
tum regimes (including the renormalized classical regime 
where, despite its name, quantum fluctuations are essen- 
tial) to the purely classical regimes, where AFMs and 
FMs exhibit identical behavior near their respective or- 
dering wavevectors. For a detailed description of the 
phase diagram and the proposed new regimes, we refer 
the reader to Fig.|l| and Table 1. 

For the antiferromagnets, both our series expansion 
calculations and the neutron data suggest that the 

regime boundaries in Fig.^ are positioned such that the 
universal T <C A scaling theory, free of any lattice cor- 
rections, is not obeyed for spin-one and larger at any nu- 
merically or experimentally accessible values of the corre- 
lation length. The data also highlights the disagreement 
in the magnitude of spin correlations even for S" = 1/2 
at all accessible temperatures with the existing RC pre- 
dictions, whereas ^(T) was found ||] to agree remarkably 
well with the RC theory. 

Our results suggest that an accurate analytical the- 
ory for S' > 1 in the experimentally relevant temper- 
ature range can be constructed by taking into account 
leading quantum corrections about the S = oo limit on 
a lattice. The important step here would be to obtain 
an analytic approximation which is valid in the classical 
crossover regime. On the other hand, any theory based 
on a purely continuum description, such as the QNLcr 
model without lattice corrections, is clearly inadequate 
for any spin larger than S = 1/2. 
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